REGIONAL, SINGLE POINT, 
AND GLOBAL BLOW-UP FOR THE FOURTH-ORDER 
POROUS MEDIUM TYPE EQUATION WITH SOURCE 



V.A. GALAKTIONOV 

Abstract. Blow-up behaviour for the fourth-order quasilinear porous medium equation 
with source, 

(0.1) u t = -{\u\ n u) xxxx + \u\ p - l u in MxM+, n > 0, p > 1, 

is studied. Countable and finite families of similarity blow-up patterns of the form 

u s (x,t) = (T-t)-^f(y), where y = x/(T - tf, = 

which blow-up as t — > T~ < oo, are described. These solutions explain key features of 
regional (for p — n+ 1), single point (p > n + 1), and global (p £ (1, n+ 1)) blow-up. The 
concepts and various variational, bifurcation, and numerical approaches for revealing the 
structure and multiplicities of such blow-up patterns are presented. 



1. Introduction: blow-up reaction-diffusion models 

1.1. On classic second-order blow-up models and higher-order diffusion. Blow- 
up phenomena as intermediate asymptotics and approximations of highly nonstationary 
processes are common and well known in various areas of mechanics and physics. The 
origin of intensive systematic studies of such nonlinear effects was gas dynamics (since 
the end of the 1930s and 1940s) supported later in the 1960s by plasma physics (wave 
collapse) and nonlinear optics (self- focusing phenomena). Nevertheless, it was reaction- 
diffusion theory that exerted the strongest influence on mathematical blow-up research 
since the 1970s. It is not an exaggeration to say that reaction-diffusion theory proposed 
basic and nowadays canonical models, which eventually led to qualitative and rigorous 
description of principles of formation of blow-up and other singularities in nonlinear PDEs. 

Finite-time blow-up singularities lie in the heart of several principal problems of PDE 
theory concerning existence, uniqueness, optimal regularity, and free-boundary propa- 
gation. The role of blow-up analysis in nonlinear PDE theory will increase when more 
complicated classes of higher-order degenerate parabolic, hyperbolic, nonlinear dispersion, 
and other equations of interest are involved in the framework of massive mathematical re- 
search and application. For such general classes of equations with typically non-potential, 
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non-divergent, and non-monotone operators (see classic monographs by Berger [2] and 
Krasnosel'skii-Zabreiko [3T] for fundamentals of nonlinear operator theory), applications 
of many known classic techniques associated with some remarkable and famous specific 
PDEs become non-aplicable, so that a principally new methodology and even philosophy 
of nonlinear PDEs via blow-up analysis should play a part. 

Reaction-diffusion models with blow-up. The second-order quasilinear reaction- 
diffusion equations from combustion theory are widely used in the mathematical literature 
and many applications. This class of parabolic PDEs of principal interest in the twenti- 
eth century includes classic models beginning with the Frank- Kamenetskii equation with 
exponential nonlinearity (solid fuel model, 1938 [31]), 

u t = Am + e n , 

u t = Au + u p , 

(1.1) 

u t = A{u n+1 ) +u p , 

u t = V ■ {\Vu\ n Vu) + u p , 

etc., where n > and p > 1 are fixed exponents, and similar equations with more general 
nonlinearities. Due to the superlinear behaviour of the source terms e u or u p as u — > +oo, 
these PDEs are known to create finite-time blow-up in the sense that a bounded solution 
ceases to exist and 

(1.2) sup^gjj u(x, t) — > +oo as t — > T~ < oo. 

If blow-up happens, the first key question is the behaviour of solutions as t — > T~ 
that reflects both mathematical and physical-mechanical essence of these phenomena. 
Such singular limits create a class of one of the most difficult asymptotic problems in 
nonlinear PDE theory. The internal structure of these blow-up singularities can be rather 
complicated even for simplest combustion models in (11.11) . It is worth mentioning that a 
full classification of the blow-up patterns is still not available for the last two quasilinear 
equations in this list. 

In combustion modelling, the blow-up phenomena (jl.2p are treated as adiabatic explo- 
sions (first results were due to Todes, 1933, in an ODE model), as an extremal type of 
instabilities in such nonlinear systems. The research in the area of blow-up combustion 
processes had been essentially intensified since the 1970s, when blow-up ideas first time 
penetrated key principles of the Laser Fusion (E. Teller's famous Report of 1972 on possi- 
bility of blow-up shockless compression of a D-T drop); see surveys in [26} p. 401] and [231 
§ 1.2]. There exists a huge amount of mathematical papers devoted to various aspects of 
blow-up and other singularity behaviour in parabolic equations and systems. We refer to 
monographs [TJ Q3H E51 EH E] , where detailed structures of blow-up are mainly studied 
and described for second-order react ion- diffusion PDEs. In [33], a systematic nonlinear 
capacity approach is developed to prove finite-time blow-up for wide classes of nonlin- 
ear evolution PDEs and systems of various orders and types. In [251 Ch. 3-5], blow-up 
structures for higher-order parabolic, hyperbolic, and nonlinear dispersion equations are 
treated by means of exact solutions on invariant subspaces for nonlinear operators. 
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First study and classification of finite families of single point blow-up patterns (p > 
7i + l) for the second-order counterpart 

(1.3) m = {u n+1 ) xx + u p (u>0) 

have been known since 1970s. Basic new ideas of blow-up and heat localization for (11.31) 
were developed by Kurdyumov and his Russian School of blow-up regimes. See key 
references, main results, and history in j4TJ Ch. 4]. 

Higher-order diffusion and structure of blow-up. Much less is known still on 
formation of blow-up singularities in parabolic PDEs with higher-order diffusion that is 
also of essential demand in modern application. Even simplest PDEs such as the extended 
Frank- Kamenetskii equation in one dimension, 

(1.4) Uf U X xxx + 6 (or Ut Uxxxxxx + 6 

and their counterparts with power source 

(1.5) u t = -u xxxx + \u\ p ^u (or u t = u xxxxxx + \u\ p ~ l u), 

reveal several principally new asymptotic blow-up phenomena demanding novel mathe- 
matical approaches; see details in [9j[TE]. Similar difficulties occur for the model equation 
from the Semenov-Rayleigh-Benard problem with the leading operator of the form 

(1-6) u t = -u xxxx + (3[{u x )\ + e u 03 >0); 

see [28]. The mathematical difficulties in understanding the ODE and PDE patterns 
increase dramatically with the order of differential diffusion operators in the equations. 
Interesting regional blow-up (see a definitions below) and oscillatory properties [10] are 
exhibited by a semilinear diffusion equation with "almost linear" logarithmic source term 

(1.7) Ut = —u xxxx + u In 4 |w|. 

There are just a few examples of higher-order parabolic react ion- diffusion PDEs with 
reasonably well- understood (but still not fully proved!) blow-up singularities. On the 
other hand, as we have mentioned, the techniques for proving global nonexistence of 
solutions are already available and cover a variety of nonlinear evolution equations of 
arbitrary orders; see [34]. Note that all the above models fll.4p - fll.7l) are semilinear and 
do not admit blow-up patterns with finite propagation (i.e., with finite interfaces) and 
hence with free boundaries that are of special interest in PDE and ODE theory. 

1.2. Fourth-order PME-type equation with source: on the model, techniques, 
and main results. In the present paper, we propose and study self-similar blow-up for 
the following fourth-order porous medium equation with source (PME— 4 with source) : 



u t = A(u) = —(\u\ n u) + \u\ p l u in R x 



where n > and p > 1. For n = 0, this gives the semilinear equation (11.51) describing single 
point blow-up for all p > 1, [9]. For n > 0, blow-up phenomena and the corresponding 
mathematics for (11. 8p are unknown and clearly become more involved. 
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We consider for (11.81) the Cauchy problem with bounded compactly supported data 
(1.9) u(x, 0) = u (x) in E. 

Unlike the well-known from the 1980s thin film equations (TFEs) (see (11.151) below), 
the PDE (11.81) contains a fully divergent, monotone, and potential, differential operator. 
Moreover, (I1.8P is a gradient system and admits strong estimates via multiplication by 
(|"u| n it)t in L 2 . The PME operator is potential in the metric of H~ 2 and is also monotone, 
so local existence of a unique continuous solution follows from classic theory of monotone 
operators; see [331 Ch. 2]. Finite propagation in the PDE (11.81) is proved by energy 
estimates via Saint- Venant's principle; see various techniques in [3] and more recent 
results in [12] and also a survey in [21] for further references. 

Global nonexistence results for PDEs like (jl.8p are known from the beginning of the 
1980s; see [17]. Concerning more recent results, let us mention that any solution of the 
Cauchy problem for 

u t = -A 2 (\u\ n u) + \u\ p in R N x R + , u(x, 0) = u (x) in R N 

(recall the difference in the source term \u\ p of the equation, which is associated with the 
nonlinear capacity method) blows up in the subcritical Fujita range 

n + l<p<p = n + l + jf, provided that f u > 

(the first condition n + 1 < p is also purely technical); see [12] for a brief simple exposition 
of the method and [31] for a systematic treatment. 

It turns out that studying and proper understanding of higher-order models such as 
(II. 8p demand completely different and more difficult mathematics, and, often, the known 
mathematical methods of classic nonlinear operator theory are not sufficient to answer 
some even basic and principal questions arisen. In other cases, tremendous efforts are 
necessary to justify their applicability for the nonlinear higher-order non-potential and 
non-variational operators involved. As happened before for the classic combustion model 
(II. 3p . complicated blow-up singularities for (jl.8p hopefully are helping to initiate new 
mathematical directions that will take a long time to be proper developed. Several open 
problems appear that are stated throughout the text. Some of them are extremely difficult 
so the author just states and discusses these without specific mathematical efforts to solve. 

Local problem: oscillatory sign-changing behaviour at finite interfaces. 
It is curious that though the principal fact on the oscillatory behaviour of solutions close 
to interfaces of the pure fourth-order PME (PME— 4) 

(1-10) «. = -(|«r«)_ (n>0) 

was rigorously established (together with existence and uniqueness of the Zel'dovich- 
Kompaneetz-Barenblatt source-type solutions) by Bernis [4 J (1988) and in Bernis-McLeod 
|6J (1991), a detailed generic structure of such oscillations for (ll.lOp was not still fully 
addressed in mathematical literature. This was one of the principle open problem of local 
theory of higher-order PMEs. Therefore, we begin our analysis with this local oscillatory 
phenomenon that is also key for some similarity blow-up profiles: 
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(0) oscillatory behaviour of solutions near interfaces (Section EJ). 

Global blow-up problems: patterns of S-, LS-, and HS-regimes. Concerning 
further evolution properties of (11. 8\\ . one of the main mathematical problems and difficul- 
ties is associated with the description and classification of blow-up patterns occurring in 
finite time. We reveal three classes of similarity blow-up for (11.81) in the ranges: 

(1) p — n + 1, regional blow-up or S-regime of blow-up in Kurdyumov's terminology 
[4"lj Ch. 4] (Section [5]): this means that the solutions blows up as t — > T™ in a bounded 
localization domain; 

(ii) p > n + 1, single point blow-up or LS-regime (Section ED: the final-time profile 
\u(x,t)\ — > oo as t — > T~ at a single point and is uniformly bounded away from it; and 

(iii) p E 1), global blow-up or HS-regime (Section [7]): |w(x, t)| — > oo as t — >■ T~ 
uniformly on any compact subset in x. 

In all the three cases, we detect the first po-branch of Fo-similarity profiles, which are 
expected to be evolutionary (structurally) stable, i.e., describe generic formation of blow- 
up singularities for the PDE I1.8| ). We present no proof (it is assumed to be very difficult) 
and no numerical evidence (a PDE modelling is then necessary that is out of our goals 
and reach here), and rely on our previous experience. 

For p G (1,71 + 1], the similarity patterns are shown to be compactly supported and 
oscillatory near interfaces, while for p > n + 1 they are not. It is also crucial that, in 
the case of regional blow-up for p = n + 1, the problem for similarity profiles becomes 
variational. Therefore, in Section [5j we detect countable families of various compactly 
supported similarity solutions by Lusternik-Schnirel'man category theory [311 § 57] and 
spherical fibering approach [38, 39] . Our rather natural and practical approach is to use 
variational families of similarity profiles for p = n + 1 as the origin of p-branches for 
p ytz n + 1 appeared at the branching point p — n + 1. In other words, the continuation 
in the parameter p makes it possible to observe and classify several blow-up patterns and 
p-branches in the non- variational cases p > n + 1 and p < n + 1. In Sections El and we 
present typical p-bifurcation (branching) diagrams. For the higher-order ODEs appeared, 
several results are essentially based on refined numerical experiments (by using standard 
codes supplied Mat Lab), which often become the only tool to check global extensibility 
properties of various branches of similarity profiles. We do not think that, in a certain 
generality, all of these conclusions can be justified completely rigorously. 

It is worth mentioning another fruitful approach to use the analogy with the linear 
bi-harmonic equation 

(1.11) u t = -u xxxx + u in R x R + , 

which is obtained from ([TUD by both limits n — > and p — > 1. A countable set of 
exponential patterns for (II. lip will be described on the basis of spectral theory developed 
in Section [2j 

1.3. On other PDEs with striking similarities of the blow-up results. The present 
research and main approaches are better understood in conjunction with the results in 
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[21], where another canonical fourth-order models of the p-Laplacian operator with source, 

(1.12) u t = ~{\u xx \ n u xx ) xx + |u| p- V n > 0, p > 1, 

is studied. It turns out amazingly that several general conclusions on blow-up behaviour 
for (11. 81) and ( 11.121) coincide, but the necessary mathematics is completely different in 
many key places (e.g., the local oscillatory behaviour and the L-S variational approach), 
since the diffusion operators works in distinct functional spaces, H~ l and L 2 respectively. 
However, even existing similarities for these two nonlinear PDEs are not that easy to de- 
tect, since the mathematics involved essentially differs and the common roots of analogous 
results are rather obscure still. We will refer and use some common results in [21] when 
necessary and appropriate. The main principles of classification of blow-up similarity 
patterns and branches are claimed to remain the same for higher-order PDEs of similar 
type, e.g., for the sixth-order PME-type equation with source (PME— 6 with source) 

(1.13) u t = {\u\ n u) xxxxxx + |u| p_ V 

Moreover, further developing this paradoxical phenomenon of analogy, the author claims 
that a number of local and similarities and geometrically equivalent set of nonlinear 
eigenf unctions are available for nonlinear degenerate PDEs such as (these awkward and 
artificial examples are chosen for the sake of the argument; note that these are not that 
far from real applications currently needed) 

(1.14) Uttt (|^| ^)xxxxxxxx ~\~ |^|^ ^ Or Uttttt ( \Uxxxxx \ ^xxxxx) xxxxx (l^l^ U) xx . 

These PDEs are even hardly classified being a mixture of parabolic (mainly this), as well 
as having some features of wave, and nonlinear dispersion types. Such equations can be 
shown to exhibit blow-up patterns with oscillatory finite interfaces, where the mathematics 
and patterns formation/classification remain similar but of course more involved. 

1.4. On other higher-order PDEs with blow-up or extinction and finite in- 
terfaces. For convenience of the Reader and for completeness of our short survey on 
blow-up in nonlinear higher-order parabolic PDEs, the interface and blow-up phenomena 
are natural and most well-known for the degenerate unstable thin film equations (TFEs) 

(1.15) u t = -(\u\ n u xxx ) x - {\u\ p - l u) xx with n > 0, p > 1. 

Equations of this form admit non-negative solutions constructed by special sufficiently 
"singular" parabolic approximations of nonlinear coefficients that often lead to free- 
boundary problems. This direction was initiated by the pioneering paper [5] and was 
continued by many researchers; we refer to [321 HE] and the references therein. Blow-up 
similarity solutions of the fourth-order TFE (I1.15P have been also well studied and un- 
derstood; see [H El [T5J 021 US] , where further references on the mathematical properties 
of the models can be found. Oscillatory solutions and countable sets of blow-up patterns 
for this TFE were described in [15] . 

Interface and extinction behaviour occur for other reaction-absorption PDEs such as 

(1.16) u t = -u xxxx - |w| p_1 w 
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in the singular parameter range p G (— 1,1), where the absorption term is not 

Lipschitz continuous at u — 0. This creates interesting evolution extinction phenomena 
and various sets of similarity patterns; see [20J and references therein. 

2. Fundamental solution and spectral properties for n = 

These questions have been discussed in [211 § 3.3], so we briefly report on those auxiliary 
properties. Consider the linear bi-harmonic equation corresponding to n = 0, 

(2.1) Ut = —u xxxx in Kx R + . 
Its fundamental solution has the form 

(2.2) b(x,t) — t~*F(y), y = x/t*, where 

(2.3) BF = -F (4) + \ yF' + \ F = in R, with / F = 1. 

The rescaled kernel F = F(\y\) is radial, has exponential decay, oscillates as \y\ — > oo, 
and [14, p. 46] 

(2.4) \F(y)\ < De~ dM4/3 in R, 

for some positive constant D and d = 3 • 2~ n / 3 . The necessary spectral properties of 
the linear operator B and the corresponding adjoint operator B* can be found in [13] 
for similar more general 2mth-order operators (see also [15, § 4]). In particular, B is 
naturally defined in a weighted space L 2 (R), with p(y) = e a \y\ A, \ where a e (0, 2d) is a 
constant, and has the discrete (point) spectrum 

(2.5) <t(B) = {A, = -|, 1 = 0,1,2,...}. 

The corresponding eigenf unctions are given by differentiating the kernel, 

(2.6) ^( 1/ ) = tJij?P)( 1/ ) | 1 = 0,1,2,.... 
The adjoint operator 

(2-7) B* = -Dl - 1 yD v 

has the same spectrum (12.51) and eigenfunctions being generalized Hermite polynomials 

(2.8) ^(y) = ^E\~=o l ^y l , 1 = 0,1,2,..., 

which form a complete and closed set in L 2 p * (R), where p*(y) = ^y. The bi-orthonormality 
of eigenfunction sets holds: 

(2.9) (i>i,i/>k) = 6 lk , 

where (•, •) denotes the standard (dual) scalar product in L 2 (R). 

3. Local asymptotic properties of solutions near interfaces 
We study the oscillatory behaviour of solutions of (jl.8p close to interfaces. 
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3.1. Local properties of travelling waves (TWs). As customary, we use TWs solu- 
tions, 

(3.1) u(x, t) = f(y), y = x- Xt, 

to describe propagation properties for (11. 8p . Then f(y) solves 

(3.2) - xf = -d/r/) (4) + i/r 1 /. 

By a local analysis near the interface, at which / = (and, as we will see, also /' = /" = 
/'" = for all n > 0), it is not difficult to see that the higher-order term l/l^" 1 / on the 
right-hand side is negligible. Therefore, near moving interfaces, we consider the simpler 
equation 

(3.3) (|/| B /T = -/ for y>0, /(0) = 0, 

which is obtained on integration once. Here we set A = —1 by scaling for propagating 
waves. It is convenient to use the natural change 

(3.4) F=\f\ n f => F'" = -\F\~^F. 

We describe oscillatory solutions of changing sign of (13.41) . with zeros concentrating at the 
given interface point y = + . Q.v. the pioneering papers by Bernis and McLeod [HE]. 

Actually, we are going to give a sharp description of the behaviour of the solutions close 
to interfaces. By the scaling invariance of (13.41) . we look for solutions of the form 

(3.5) F(y) = y^(p{s), s = lny, where fj, = 3 ^ n ^ > 3 for n > 0, 

where <p(s) is the oscillatory component. Substituting ( 13.51) into ( 13.4ft yields the ODE 

n 

(3.6) P 3 (<p) = -\<p\ ^<f. 

Linear differential operators P^ are given by the recursion (see [2"5j p. 140]) 
P k +i(<p) = P' k {y) + (/i - k)P k ((p), k>0; P Q ((p) = (p, so that 
P^cp) =if' + fop, P 2 (<p) = if" + (2/i - l)<p' + M/i - l)ip, 
P 3 {(p) = ip'" + 3(/i - l)ip" + (3fi 2 - 6// + 2)(p' + /i(/i - l)(/i - 2)(p, etc. 

In view of (I3.5p . we look for uniformly bounded global solutions <p(s) that are defined 
for all s — 1, i.e., can be extended to the interface at y — + . The best candidates 
for such global orbits of (13.61) are periodic solutions <p*(s) that are defined for all s£l 
Indeed, they can describe suitable (and, possibly, generic) connections with the interface 
at s — — oo. Existence of such a periodic solution v ? *( s ) of (13.61) can be achieved by 
shooting arguments; see details in [16, § 7], where further references concerning periodic 
orbits of higher-order ODEs are given. Uniqueness of (p* remains an open mathematical 
problem, though was always confirmed numerically. 

Figure [1] shows fast convergence to such a unique (numerically) stable periodic solution 
of (13. 6p for various n > 0. These numerics are obtained by using the ODE solver ode45 
MatLab with enhanced tolerances and a regularization in the singular term up to 10~ 12 ; 
see ( I5.17P and more comments in Section [5j where more advanced numerical techniques 



(a) n = 0.75 



(b) n = 1 




for solving boundary-value problems are employed. Different curves therein correspond 
to different Cauchy data (p(0), <£>'(0), ^"(0) prescribed at s — 0. For n smaller than |, 
the oscillatory component gets extremely small, so an extra scaling is necessary, which is 
explained in [16, § 7.3]. A more accurate passage to the limit n — > in (13.61) is done there 
in § 7.6 and in Appendix B. In (d), we also present the periodic solution for n = +oo 
where (13 .6p takes a simpler form, 

p 3(p) = -signy?. 

Finally, given the periodic of (I3.6p . as a natural way to approach the interface 

point y = according to (13.51) . we have that the ODE (13.31) and, asymptotically, (13. 2\f , 
generate at the singularity set {/ = 0} 

(3.7) a 2D local asymptotic bundle with parameters yo and phase shift in s i— > s + so- 
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This two-dimensional bundle will be shown to be sufficient for matching with, typically, 
two boundary conditions that generate blow-up patterns. 

3.2. Non-oscillatory case A > 0: ID bundle of non-oscillatory asymptotics. For 

A = +1, we have the opposite sign in the ODE 

(3.8) Ps(<p) = \<p\ 
which admits two equilibria 

n + 1 

(3.9) <p± = ±[ii(ji-l)(ji-2)] » . 

Let us check the dimension of their stable manifolds as s — > +oo and, most importantly, 
s — ► — oo, which corresponds to approaching the interface. Note that, for (13. 6p the trivial 
equilibrium (po — is evidently unstable in both directions with empty stable manifolds 
(in view of the non-Lipschitz nonlinearity at (p — on the right-hand side). 

Thus, by standard linearization, it follows that both are stable as s — > +oo: setting 
if = <£>+ + Y yields the linear ODE 

(3.10) PY = [p 3 --i T ^-l)( /U -2)/]Y = 0. 

The characteristic equation by plugging K = e As into (13 . 10[) takes the form 

(3.11) p(A) = A 3 + 3(/i- 1)A+ (3/i 2 - 6/x + 2)A + 3(// - 1)0* -2) = 0, where 

p(0) = 3(> - 1)0 - 2) > 0, j/(A) = at A± = f [ - 3{pt - 1) ± y/3] < 0. 
It follows that all characteristic values of P satisfy 

(3.12) ReA fe <0, jfe = l,2,3. 

Figure [2](a) also confirms that as s —>■ +oo the equilibria ( 13.91) are stable. In (b), which 
gives the enlarged image of the behaviour from (a) close to if = 0, we observe a changing 
sign orbit, which is not a periodic one. Therefore, this cannot be extended as a bounded 
solution up to the interface at s = — oo. There are many other similar ODEs, where such 
behaviour between two equilibria is periodic; cf. [251 P- 143]. 

In a whole, theses results confirm that for A > 0, the TWs are not oscillatory at 
interfaces, so that, formally, this backward propagation can be performed by positive 
solutions. On the other hand, (13.121) establishes that the stable manifold of equilibria <p± 
as s — » — oo (i.e., towards the interface) is empty, so that unlike (13.71) . for A = +1, 

(3.13) 3 a ID local asymptotic bundle with parameter y . 

One can see from (13.51) that these positive asymptotics are given by (for fixed A = +1 
which has been scaled out as a non-essential parameter) 

(3.14) f(y) = O-^t^O- 1)0 -2)]-*(l + o(l)) as y - y+ . 

This ID bundle is not enough for construction of typical global connections via two 
boundary conditions, so that the backward propagation is either not possible at all for 
almost all (a. a.) initial data in the Cauchy problem for (11.101) . or it is not performed by 
TWs. 
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s s 

(a) Stability of equilibria (b) Behaviour in between 



FIGURE 2. Non-oscillatory behaviour for the ODE (|3,8p for n = 5; stability of 
equilibria (|3.9p (a), and enlarged non-periodic behaviour in between, (b). 



In other words, this shows that, for the PME— 4 (ll.lOp . the forward propagation of 
interfaces is common. Concerning the backward one, though existing for specific initial 
data, it is not plausible in general. Recall that for the standard PME 

(3-15) ut = (\u\ n u) xx , 

the backward propagation of non-negative solutions is completely prohibited by the Max- 
imum Principle (by the straightforward comparison from below with slowly moving small 
TWs). For higher-order equations, these barrier techniques via comparison fail, but any- 
way, via the reduction of the bundle dimension in (I3.13p . we justify a similar phenomenon, 
which is now true not for all, but for a. a. initial data. It is worth mentioning that the 
PME (I3.15P also admits oscillatory solutions near interfaces (see references and comments 
in [T9t p. 30]). By Sturm's First Theorem on non- increase of the number of zeros, such 
behaviour is not generic for second-order parabolic PDEs in the sense that such solutions 
cannot appear from data with any finite number of sign changes. 

4. Blow-up similarity solutions: problem setting and preliminaries 

4.1. Posing the similarity problem. The parabolic PDE (11.81) admits the following 
similarity solutions with blow-up as t — > T~: 

(4.1) u s (x,t) = (T-t)-^f(y), y = x/(T-ty, with (5 = 
The rescaled profile f(y) solves a quasilinear fourth-order ODE 

(4.2) A (/) = -(|/r/) (4) -/W' -^ T / + l/r 1 / = in E. 
Equation (14. 2 p has the constant equilibria 

(4.3) ± /,(p) = ±(p - ±oo as p->l+. 
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For even solutions f(y), symmetry conditions at the origin y = are posed, 

(4.4) f(0)=0 and /'"(0) = (if /(0) ^ 0), 
while for odd ones we impose the anti-symmetry ones 

(4.5) /(0) = and f"(0) = 0. 

A natural setting for the Cauchy problem assumes that, for p e (1, n + 1], 

(4.6) /(?/) is sufficiently smooth and compactly supported. 

The actual regularity of f(y) close to interfaces was determined in the previous section by 
the local asymptotic analysis. For p > n + 1, the solutions are not compactly supported. 

The ODE ( 14.21) is a difficult fourth-order nonlinear equation, so that, in general, exis- 
tence and multiplicity results for the above boundary-value problems reduce to compli- 
cated, at least, 2D shooting problems (recall (13.71) and two boundary conditions). This 
principally differs (14.21) from its second-order counterparts associated with the combustion 
model (11.31) . where the most of the results admit either a clear phase-plane interpretation 
or reduce to ID shooting in the presence of the Maximum Principle; see typical results 
and references in [4"Tl Ch. 4]. Bearing in mind that, for the PME— 6 with source (I1.13p . 
the parameter space of shooting will be 3D, we cannot rely on these geometric ideas, and 
have to use a principally different approach for reliable detecting multiple blow-up pat- 
terns. This will be variational approach to the case p = n + 1 plus p-branching approach 
for p > n + 1 and p < n + 1. 

4.2. Blow-up self-similar profiles: preliminaries. We next study global compactly 
supported solutions of the ODE (14. 2p . For p < n + 1, the local interface analysis from 
Section [3] applies to (14.21) . Indeed, close to the interface point y = yo > of the similarity 
profile f(y), the ODE (14. 2 p contains the same leading terms as in (13.31) and other linear 
two are negligible as y — > y$. 

It is key that, taking into account the local result (13. 7p and bearing in mind the two 
boundary conditions in (I4.4p or (14. 5p . we may expect that 

(4.7) there exists not more than a countable set of solutions. 

This speculations assume a certain "analyticity" hypothesis concerning the dependence 
on parameters in the degenerate ODE (14.21) . which is plausible but not easy to prove. 
Actually, this means that, relative to the parameter p > 1, we can expect at most a 
countable set of p-branches of solutions. To begin with, this is true for the linear case 
n = and p = 1 . 

4.3. Countable set of similarity solutions for n = 0, p = 1, and homotopic 
connection. These questions are discussed in detail in [2TI § 3.4], so we briefly comment 
that the linear PDE (11.111) has a countable set of patterns given by eigenfunctions (12. 6p : 

(4.8) ui{x,t) =e- t r 1 i i i; l (^ I ), 1 = 0,1,2,... 

Our intention is to see how the blow-up similarity patterns (14. ip can be deformed as n — > 
and p — > 1 to those in (14.81) . Then (14.81) suggests that, for n > 0, there exists a countable 
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number of "branches" {fi(y, n,p)}, which "bifurcate" from the point {n = 0, p = 1}. It is 
possible to observe the space-time structures in (14.81) in nonlinear blow-up analysis, where 
similar countable sets of patterns will be shown to exist; see |29j. We then say that the 
above two (linear for n = 0, p = 1 and nonlinear for n > 0, p > 1) asymptotic problems 
admit a continuous homotopic connection as n — > 0, p — > 1. See [21] for extra details. 

5. Regional blow-up similarity profiles for p = n + 1 

As in [211 § 4], the case p = n + 1 is simpler since (3 = 0, so (14.11) is a solution in 
separable variables 

u{x,t) = (T-t)-*f(x). 
Then / = f(y) (we continue to use y as the independent spatial variable for application 
to p 7^ n + 1) solves an autonomous fourth-order ODE of the form 

(5.1) A(/) = -(|/|"/)( 4 )-i/+|/r/ = in R. 
We again use the same change as in (13.41) . 

(5.2) |/| n / = rT^F =^ F (4) = F- \F\~^F in R. 

This equation is much simpler than that in [211, anc ^ was studied in [23]. We briefly 
present the main conclusions that will play an important role for the cases p > n + 1 and 
p < n + 1 later on. 

5.1. Variational setting. The operators in (15.21) are potential, so the problem admits a 
variational setting, so the solutions can be obtained as critical points of a C 1 functional 
of the form 

(5.3) E(F) = -y(F"ydy + yF*dy-lJ\F\»dy, where „ = ^ e (l,2). 

We are interested in critical points in W / 2 2 (R)nL 2 (R)nL I/ (R). Especially, we are interested 
in localized compactly supported solutions, so we choose a sufficiently large interval B R = 
(—R,R) and consider the variational problem for (15.31) in W^iB^), where we assume 
Dirichlet boundary conditions at the end points 8Br = {±R}. By Sobolev embedding 
theorem, W^q^B^) C L 2 (Br) and in L p+1 (5^) compactly for any p > 1. Continuity of 
any bounded solution F(y) is guaranteed by Sobolev embedding if 2 (R) C C(R). We also 
need the following: 

Proposition 5.1. Let F be a continuous weak solution of the equation (15.^) such that 

(5.4) F(y) ^0 as y -> oo. 
Then F is compactly supported in R. 

Proof. Consider the corresponding parabolic equation with the same elliptic operator, 

(5.5) w t = —w xxxx + w — \w\ p ~ 1 w in R + x R, where p = ^-j- G (0, 1), 
with initial data F(y). Setting w = e l w yields the equation 

(5.6) w t = -w xxxx - e »+i \wf l w, 
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where the operator is monotone in L 2 (M). Therefore, the CP with initial data F has a 
unique weak solution, [321 Ch. 2]. Thus, (I5.5P has the unique solution w(y,t) = F(y). 
In the presence of the singular absorption — |w| p-1 -u, with p < 1, there occurs the phe- 
nomenon of instantaneous compactification or shrinking of the support of the solution for 
any data satisfying (15.41) (or even for more general data in L p -spaces). Such phenomena 
for quasilinear absorption-diffusion equations for p < 1 have been known since the 1970s. 
By energy estimates, similar results were proved for a number of quasilinear higher-order 
parabolic equations with non-Lipschitz absorption terms, [3j H2]. By the instantaneous 
compactification, the multiplier e "+ 1 in the absorption in (I5.6P changes nothing. □ 

Thus, to revealing compactly supported patterns F(y), we have to pose the problem in 
bounded sufficiently large intervals. 

5.2. L— S theory and direct application of fibering method. Unlike J2Tj § 4], this 
application is more standard. Namely, we apply classic Lusternik-Schnirel'man (L-S) 
theory of calculus of variations [3TT § 57] in the form of the fibering method [381 [39] . 
Then the number of critical points of the functional (15.31) depends on the category (or 
Krasnosel'skii's genus [311 § 57]) of the functional subset, on which fibering is taking 
place. The critical points of E(F) are obtained by the spherical fibering in the form 

(5.7) F = r(v)v (r > 0), 
where r(v) is a scalar functional, and v belongs to the subset 

(5.8) H = {veWl (B R ): H (v) = -J(v") 2 dy + Jv 2 dy = l}. 
The new functional 

(5.9) H(r,v) = \r 2 - \r v § \v\ v dy 
has the absolute minimum point, where 

(5.10) H' r = r - r^ 1 J \v\ u dy = r (v) = (J \v\ v dy) ^ , 
at which H(r (v),v) = —^ L r 2 (v). Introducing 

(5.11) H(v)= [-£. H (r (v),v)]*? =J\v\ v dy, 

yields an even, non-negative, convex, and uniformly differentiable functional, to which 
L-S theory applies, [2U § 57]; see also [TTJ p. 353]. Following [55] . searching for critical 
points of H in 7i , one needs to estimate the category p of the set Ho. The details on 
this notation and basic results can be found in Berger [2, p. 378]. Notice that the Morse 
index q of the quadratic form Q in Theorem 6.7.9 therein is the dimension of the space 
where the corresponding form is negatively definite. This includes all the multiplicities of 
eigenf unctions involved in the corresponding subspace. 

For application, it is convenient to recall that utilizing Berger's version [21 p. 368] of 
this minimax analysis of L-S category theory [311 p. 387], the critical values {c^} and the 
corresponding critical points {t^} are given by 

(5.12) c k =mi T£Mk sup^ H(v), 
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where T C 7io are closed sets, and M. k denotes the set of all subsets of the form BS^ 1 C 
TCo, where S k ~ l is a suitable sufficiently smooth (k — l)-dimensional manifold (say, sphere) 
in 7io and B is an odd continuous map. Then each member of M. k is of genus at least k 
(available in 7io)- It is also important to remind that the definition of genus [3U p. 385] 
assumes that p(J-") = 1, if no component of T U J 7 *, where J 7 * = {v : —v G J 7 }, is 
the reflection of T relative to 0, contains a pair of antipodal points v and v* = —v. 
Furthermore, p(.F) = n if each compact subset of T can be covered by, minimum, n sets 
of genus one. 

According to (15.121) . c\ < c 2 < ... < Q , where l = Iq(R) is the category of TCq (see an 
estimate below) satisfying 

(5.13) Iq(R) +oo as R — ► oo. 

Roughly speaking, since the dimension of the sets JF involved in the construction of M. k 
increases with k, this guarantees that the critical points delivering critical values (15.121) 
are all different. 

It follows from (15.81) that the category l = p(Ho) of the set Ho is equal to the number 
(with multiplicities) of the eigenvalues X k > — 1 of the linear bi-harmonic operator, 

(5.14) -w^ = X k ^ } 1>eWi fi (B R ). 
Since the dependence of the spectrum on R is, obviously, 

(5.15) X k (R) = iT 4 A fc (l), k = 0,1,2,..., 

the category p(Tlo) can be arbitrarily large for R ^> 1, (15.131) holds, and we obtain: 

Proposition 5.2. The ODE problem (\5.2$ has at least a countable set of different solu- 
tions denoted by {Fi, I > 0}, each one obtained as a critical point of the functional (\5.3§ 
in W^ (Bf() with sufficiently large R > 0. 

5.3. First basic pattern and local structure of zeros. We next present numerical 
results concerning existence and multiplicity of solutions for equation (15.21) . In Figure [31 
we show the first basic pattern for (15.21) called the Fo(y) for various n e [0.1, 100]. Note 
that (I5.2p admits a natural passage to the limit n — > +oo that gives the ODE 

F - 1 for F > 0, 



(5.16) F {4) = F-signF = | 



F + 1 for F < 0. 



A unique oscillatory solution of (15.161) can be treated by an algebraic approach; cf. 
§ 7.4]. The solutions for n = 1000 or n = +oo practically do not differ from the last 
profile for n = 100 in Figure [3J 

These profiles are constructed by MatLab by using a natural regularization in the sin- 
gular term in (15. 2p . 

(5.17) F^ = F - (e 2 + F 2 ) ^) F in R. 

Here, the regularization parameter e and both absolute and relative tolerances in the 
bvp4c solver have been enhanced and took the values 

(5.18) e = 10~ 10 and Tols = 10~ 10 . 

15 



This allows us also to reveal the refined local structure of multiple zeros at the interfaces. 
Figure H] for n = 1 shows how the zero structure repeats itself in a self-similar manner 
from one zero to another in the usual linear scale. 

In Figure [5j we present the results for n = 1 that show the oscillatory structure in the 
log-scale such as (13.51) with y t— > y — y and 

(5.19) fi = ^| n=1 = 8 =► In \F(y)\ = 8\n(y -y) + \n \^(\n(y - y))\ + ... 

(see also (15.231) below) . This figure shows an "e-dynamic" formation of at least six-seven 
nonlinear zeros, when we decrease e = Tols from 10~ 3 (just two first nonlinear zeros; the 
rest correspond to the linear ones in (I5.20p ) to 10~ 10 (6-7 zeros are nonlinear). Observe a 
"concave" shape of the graph for the last e = 1CT 10 , which is consistent with the log-shape 
in (I5.19P (we claim that even a trace of the multiplier 8 can be distinguished by 6 zeros). 

Further decrease of e and Tols leads to quick divergence of bvp4c solver, which has the 
limit minimal Tols 10~ 13 , which is not that helpful in comparison with 10~ 10 achieved in 
Figure [5j The right-hand interface is then estimated as (it is quite a challenge to detect 
numerically the free-boundary point more accurately) 

y « 12. 

The "nonlinear area" ends at y ~ 11 (due to the requires accuracy and £-regularization) , 
and next we observe the "linearized area" where (I5.17P implies an exponential behaviour 
for y ^> 1 governed by the linearized ODE 

— 71 

(5.20) F (4) = -e'^F + ... F{y) ~ e -ci^ 4( " +1) cos (c 2 e~^^y + c 3 ) , 

where Ck are constants and c\ > 0. From y ~ 12.5 (F ~ 10 -12 ), the numerics are not 
reliable at all. 

5.4. The first basic pattern F (y): on the behaviour as n — » 0. In order to under- 
stand the continuous connection with the linear PDE (11. lip and other homotopy results, 
one needs to pass to the limit n — > + in the ODE (I5.2p . Figure [6] shows the behaviour of 
the first basic profile F (y) for small n = 10 _1 , 10~ 2 , ... , 10~ 7 (a), while (b) explains the 
behaviour of the maximum value -Fo(O) and around. 

It follows from equation f)5.2p that, for n k + , close to the origin and uniformly on 
bounded intervals in y, the smooth solution F (y) approaches an even function satisfying 
the linear homogeneous equation 

(5.21) F (4) =0 F (y) w C (n) - A (n)y 2 , 

so that F (y) essentially depends on small n. Here both unknown coefficients Cq(ti) — > 
C*o(0) and A (n) — > are slightly oscillatory as n — > (in view of the non-monotone 
behaviour in (15.231) near the interface), with 



(5.22) 



C (0) w 1.435. 
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Basic patter F Q (y) for various n 



n=100 






n=5 



























n=1 



n=0.5 



FIGURE 3. The first solution Fq(u) of (|5.2p for various n. 



_ 3 m=2, n=1: basic pattern F (y), oscillations enlarged 




_ 4 m=2, n=1 : basic pattern F (y), oscillations enlarged 




(a) scale 1(T 3 (b) scale 1CT 4 

FIGURE 4. Enlarged zero structure of the profile Fq(jj) for n = 1 in the linear scale. 



The behaviour f 1 5 . 2 1 1) is then to be matched with the asymptotics near the interface at 
some y = y$ (n) ^> 1 (cf. (13.51) ) that is governed by the non-Lipschitz term in ( 15.21) . 

IA\ . n , , 4(n + l) 

(5.23) F^ = -\F\ + =^ F(y) = (y -y) » <p m ( ln(y - y)) + - , 

where is the corresponding periodic solution of the ODE for the oscillatory component. 
A dimensional analysis implies that the interface position can be estimated as follows: 

(5.24) yo(n) ~ as n — * + ; 
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FIGURE 5. Behaviour of -Fo(y) for re = 1 in the log-scale near the interface: 
formation of five nonlinear zeros as e = Tols decrease from 10 -3 to 1CP 7 . 



n-deformation of F Q (y) from n=1 (dashed line) to n=10 6 Enlarged: F Q (y) from n=1 (dashed line) to n=10' 




-200 -150 -100 -50 50 100 150 200 -10 -5 5 10 



y y 

(a) profiles Fo(y) (b) enlarged for y w 

FIGURE 6. Deformation of the similarity solutions F${y) of the ODE (|5.2p for 
small n > 0. 

see [ini § 9.6] for details. First rough matching of (15.21 1) and (15.231) . (15.241) yields 

A (n) ~ «2. 

A more accurate matching is difficult and typically assumes appearing logarithmic, Inn, 
and other terms including oscillatory ones; see [161 § 7]. 



is 



5.5. Basic countable family. Figure [7J showg[j] the basic family denoted by {Fi, I = 
0,1,2,...} of solutions of (15. 2p for n = 1. This is associated with the application of 
L-S and fibering theory, [23]. Each profile Fi(y) has / + 1 "dominant" extrema and I 
"transversal" zeros; see [23], § 5] and [221 § 4] for further details. We claim that 

all the internal zeros of Fi(y) are transversal 

(excluding the oscillatory end points of the support), and this is a fundamental conclusion 
of our analysis, which deserves more rigorous treatment. In other words, each profile Fi 
is approximately obtained by a simple "interaction" (gluing together) of I + 1 copies of 
the first pattern ±F Q taking with necessary signs; see further comments below. 

Actually, if we forget for a moment about the complicated oscillatory structure of solu- 
tions near interfaces, where an infinite number of extrema and zeros occur, the dominant 
geometry of profiles in Figure [7J approximately obeys Sturm's classic zero set property, 
which is true rigorously for the second-order ODE only, 

(5.25) F" = -F+\F\ n+1 F in R. 

For (15.251) . the basic family {Fi} is constructed by direct gluing together simple patterns 
±F Q given explicitly; see [25j p. 168]. Therefore, each Fi consists of I + 1 patterns (with 
signs ±i*b), so that Sturm's property is clearly true. 

5.6. Countable family of {F , F }-gluing and extensions. This procedure is similar 
to that in [2Tj § 4.7], but we again recall that these are completely different variational 
problems. Further patterns to be introduced do not exhibit as clear a "dominated" 
Sturm property and are associated with a double fibering technique where both Cartesian 
and spherical representations are involved; see details [231 § 6]. Let us present basic 
explanations. 

The nonlinear interaction of the two first patterns Fo(y) leads to a new family of profiles. 
In Figure [8] for n — 1, we show the first profiles from this family denoted by {F +2t k,+2}, 
where in each function F +2 ,k,+2 the multiindex a = {+2, k, +2} means, from left to right, 
+2 intersections with the equilibrium +1, then next k intersections with zero, and final 
+2 stands again for 2 intersections with +1. Later on, we will use such a multiindex 
notation to classify other patterns obtained. 

Since F Q (y) is infinitely oscillatory at the interfaces, the family {F +2i a:,+2} is expected 
to be countable, so such functions exist for any even k = 0,2,4, ... . Therefore, k = +oo 
leads to the non-interacting pair with no overlapping of supports, 

(5.26) F (y + yo) + F (y - y ), where supp F (y) = [-y , y }. 

There exist various triple {F ,F ,F } and any multiple interactions {F ,...,F Q } of arbi- 
trarily large k single profiles, with a variety of distributions of zeros between any pair of 
neighbours. 

1 Amazingly, these profiles look almost identical to those for the PDE ()1.12|) [21| § 4], corresponding 
to a completely different diffusion operator and much harder and distinct L-S theory. These similarities 
underline some quite obscure and involved common evolution properties of blow-up for different PDEs. 
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y 

FIGURE 8. First patterns from the family {-F+2^+2} of the {Fq, Fo}-gluing; n = 1. 



m=2, n=1: family {F_ 2k+2 } 




y 

FIGURE 9. First four patterns from the family {F_2,fe,+2} of the {— Fo,Fo}- 
gluing; n = 1. 

5.7. Countable family of {— Fq, Fo}-gluing and extensions. Consider next the in- 
teraction of —F (y) with F (y). In Figure [91 for n — 1, we show the first profiles from this 
family denoted by {F_ 2 ,k,+2] , where, for the multiindex a = {—2, k, +2}, the first number 
—2 means 2 intersections with the equilibrium —1, etc. The enlarged zero structure shows 
that the first two profiles belong to the same class -F-2,i,2 ; i-e., both have a single zero 
for y w 0. The last solution shown is F_ 2)5i+2 . This family {F_ 2 ,fc,+2} is expected to be 
countable, with profiles existing for any odd k — 1, 3, 5, where the pair for k = +oo is 
non-interacting, — F (y +y ) + F (y — y ). There can be constructed families of an arbitrary 
number of interactions such as {±Fo, ±Fq, ±Fo} consisting of any k > 2 members. 
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m=2, n=1 : Unstable periodic solution about 1 : F(0)=1 .5, F"(0)=-0.3787329255... 
2 
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0.2 


5 10 15 20 25 30 

y 

FIGURE 10. An example of a periodic solution of the ODE (|5.2p for n = 1. 

5.8. Periodic solutions in R. We are going to introduce new types of patterns, which 
exhibit a different geometric shape. We then need first to describe non-compactly sup- 
ported periodic solutions in R. As a variational problem, equation (15.21) admits an infinite 
number of periodic solutions; see e.g. [3H Ch. 8]. Figure [TD] for n — 1 presents a special 
unstable (as y — ► +oo) periodic solution obtained by shooting from the origin with condi- 
tions F(0) = 1.5, F'(0) = F"'(0) = 0, and F"(0) = -0.3787329255... . The periodic orbit 
F*(y) with the value 

(5.27) F*(0) « 1.535... 

will be shown to play a key role in the construction of other families of compactly sup- 
ported patterns. Namely, this variety of solutions of (I5.2p having oscillations about equi- 
libria ±1 are close to ±F*(y) there. 

5.9. Family {F + 2k}- The patterns F + 2k for k > 1 have 2k intersections with the single 
equilibrium +1 only and have a clear "almost" periodic structure of oscillations about it; 
see Figure ITTT a). The number of intersections with F = +1 denoted by +2k is an extra 
characterization of the Strum index to such a pattern. In this notation, F + 2 = Fq. 

5.10. Complicated and chaotic patterns. On the basis of our previous experience of 
dealing with various patterns, we will classify other solutions (possibly, a class of patterns) 
by introducing multiindices of any length 

(5.28) a = {±01, a 2 , ±cr 3 , 04, ±07}. 

Figure [TlT b) shows several profiles from the family with the index a = {+k, I, —m, I, +k}. 
In Figure [T2J we show further two different patterns, while in Figure [T31 a single most 
complicated pattern is presented, for which 




(5.29) 



a = {-8, 1, +4, 1, -10, 1, +8, 1, 3, -2, 2, -8, 2, 2, -2}. 
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(a) F +2k (y) 

FIGURE 11. Two families of solutions of (|5.2p for n 

F+k,l,—m,l,+k 

(b). 



(b) f+k,i, 

= 1; F+2k(y) (a) and 



m=2, n=1 : pattern F 



+2,2,+4,2,+2,1,-4 



m=2, n=1 : pattern F 



+6,3,-4,2,-6 




10 20 




(a) cr = {+2, 2, +4, 2, +2, 1, -4} (b) a = {+6, 3, -4, 2, -6} 

FIGURE 12. Two patterns for (|oT2j) for n = 1. 



All computations are performed for n = 1 as usual, and numerics were well converged with 
sufficient accuracy and regularization at least ~ 10~ 4 — 10~ 5 and better. This shows that 
the multiindex (I5.28P can be arbitrary, i.e., can take any finite part of any non-periodic 
fraction. Though we do not insist that, for a given a, the profile F a (y) is unique, we have 
seen that the uniqueness fails. Note that the homotopy approach [5DT HI] does not apply 
to ODEs such as (15.21) with infinite oscillatory properties. 

Actually, this means chaotic features of the whole family of solutions {F a }. These 
chaotic types of behaviour are known for other fourth-order ODEs with coercive operators, 
[361 P- 198]. 
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m=2, n=1: pattern F^^^^,,^ 



F(y)i 




FIGURE 13. A complicated pattern F a (y) for ([P]) for n = 1. 

6. Single point blow-up for p > n + 1 

We now consider the ODE (14.21) in the case p > n + 1, which, in view of the spatial 
rescaled variable y in (14.11) . corresponds to single point blow-up. It is crucial that (14.21) 
for is not variational. Therefore, the solutions of ( 14. 21) can be traced out by a 

complicated shooting and matching procedures, which are still not completely justified. 
For practical reasons, we will use a continuation in parameters approach, which allows us 
to predict solutions by using those in the variational case p — n + 1. Recall that for such 
ODEs, using a standard inverse function theorem is not straightforward at all since the 
differential operator in (14. 2p is degenerate and singular. Nevertheless, we will suggest to 
use Schauder's fixed point theorem and arrive at some convincing conclusions concerning 
the solvability and the multiplicity of solutions (the so-called p-branches of solutions). 

6.1. Asymptotics at infinity and first numerical results. Recall that, according 
to our local analysis close to interfaces, such a behaviour for the TWs given by (I3.3|) . 
(13 .4p asymptotically coincides with that for the blow-up similarity solutions; cf. ODEs 
(I5.2p and (15.231) . We can continue to use this convenient analogy in the present case, 
since the ODE (14. 2 p provides us with similar asymptotics at the interfaces. Obviously, 
the behaviour at finite interface yo > now corresponds to A = +1, and therefore the 
one- dimensional bundle (13.131) is not sufficient to shoot two boundary conditions in ( 14.41) 
or (14. 5p . This explains why we need another asymptotic expansion as y — > +oo. Thus, 
compactly supported similarity profiles are very unlikely (though may exist for some very 
special parameters values involved). 

Thus, unlike the previous case of regional blow-up for p — n + 1, in the present case, in 
order to get a guaranteed successful shooting, we must use another full two-dimensional 
bundle of non-compactly supported solutions of (14.21) with the following behaviour: 

(6.1) f(y) ~ W + 
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as y — > +oo. 



n=1 , p=2, 2.25, 2.50, 2.75, 3, 4, 5, 6, 7 
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FIGURE 14. Single point blow-up patterns F (y) of (|4.2p for n = 1 and various 
P > 2. 

Here Co 7^ and C\ G R are arbitrary constants and 

7 = rzTT<0, t/ = ,r 4(p rSM >° ; b o = i (4 T C(7 n ) 1 > 0. 

' p— (n+l) ' 3[p— (n+1)] ' u v Vn+l u / 

Roughly speaking, the first bracket in (16. 1 p represents an "analytic" part of the expansion 
(e.g., for integer p, in terms of multiples of y 1 it can be represented as an analytic series; 
the proof of convergence is difficult), while the second braces gives the essentially "non- 
analytic" part. Such a structure in (16 .ip is typical for saddle-node equilibria, [37J P- 311]. 
Of course, for the fourth-order ODE (14.21) . this expansion does not admit a simple phase- 
plane interpretation (though the algebraic origin of the expansion is clear). Justification of 
(16. ip needs technical applications of fixed point theorems in weighted spaces of continuous 
functions defined for y ^> 1; see a typical example in [4Tj p. 29] and related references. 

One can see passing to the limit t — > T~ in (14.11) that the asymptotic behaviour (16.11) 
gives the following final-time profile of this single point blow-up for even profiles / = 

(6.2) u s (x,T~) = Co\x\~v-^+v < 00 for all x ^ 0. 
Returning to the asymptotic expansion, we conclude that (16.11) represents 

(6.3) a 2D asymptotic bundle. 

Hence, the bundle (16.11) is well suitable for matching with also two symmetry conditions 
at the origin (14.41) . so we expect not more than a countable set of solutions. For first 
patterns, we keep the same notation as in Section [S]) for p — n+ 1. Moreover, one can 
expect that these profiles can be continuously deformed to each other asn-> + . 

In Figure we present the first pattern F (y) for n — 1 with p = 2 (the dotted line 
for comparison), 2.25, 2.5, 2.75, 3, 4, 5, 6, 7. This shows that, for larger p, the profiles get 
the positive asymptotic behaviour (16. ip with Co > 0, and become strictly positive in R. 
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Profiles F_0 and F_1 , and three of \{-F_0,F_0}-gluing; n=1 , p=2.25 
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FIGURE 15. Single point blow-up patterns F (y) and F\{y) of (|4.2p and three 
profiles from the family {-F-2,2fc+l, +2} of the {— Fq, Fo}-gluing; n = 1, p = 2.25. 

Figure [151 shows the first F (y) and the second Fi(y) (dipole-like) patterns for n = 1 and 
p = 2.25. In addition, we show therein three profiles from the corresponding (possibly, 
finite for p > n + 1) family of {— F , Foj-gluing. This family can be viewed as the 
continuous p-extension of the family shown in Figure M in the variational case p — n+ 1. 
Observe that we have two profiles F_ 2 ,3,+2(y) , and in the last one F_2,5,+2(y) the structure 
of zeros close to y = was not identified completely clear. Numerically, this is very 
difficult and convergence is very slow with the maximal number 20 000 of mesh points 
used in bvp4c solver with Tols ~ 3 • 10 -3 achieved. 

These numerics also confirm that the profiles have non-oscillatory asymptotic behaviour 
as in (16. ip . though these can change sign a few (couple of) times, thus inheriting this finite 
oscillation property from the oscillatory one of the type (13.51) for p = n+ 1 = 2 (sufficiently 
close to the current p = 2.25). 

6.2. On branching of solutions from variational critical points: degree theory 
and Schauder's theorem. We again reduce (14.21) to the semilinear problem, 

(6.4) F = \f\ n f => -F^ - (3{l- a)\F\- a F'y - ^\F\~ a F +\F\P {1 ~ a ^ l F = 

where a = For p = n + 1, i.e., for (3 = 0, this ODE has been studied in Section [5j 
Setting e = p — (n + 1) and writing (16.41) as 

(6.5) G(F) = -F^ - \F\- a F + IF^-^F = e \F\~ a F'y, 

we have that the non-autonomous term on the right-hand side becomes an asymptotically 
small perturbation for \s\ 1 (in the "weak" integral sense for the equivalent integral 
equation, so the singularity \F\~ a F' is supposed to be eliminated via integration by parts; 
see (16.61) ). On the other hand, the operator on the left is variational and hence the non- 
perturbed problem for e = admits families of solutions described in Section [5] (see more 
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details in [23]). Classic perturbation and branching theory [31] suggests that, under 
natural hypotheses, the variational problem for e = can generate a countable family 
of p-branches, which can be extended for some sufficiently small \e\ > 0. The analysis 
of bifurcation, branching, and continuous extensions is then performed for the equivalent 
to (16.51) integral equations with Hammerstein compact operators; see typical examples 
in [9j [22l 29J, where similar perturbation problems for blow-up and global patterns were 
investigated. Note that these p-branches of solutions are fully extensible and can end up 
either at a singularity point or at another bifurcation values. 
Thus, for convenience, we write (16.51) as follows: 

(B* - I)F = h(F, e) = -l y F' + e^ \F\~ a F'y + ^ \F\-*F 
_\ F \p(i-a)-i F =^ F — A(F, e) = (B* — iy l h(F, e), 

where (B* — I)^ 1 is a compact linear operator in L^»(IR), [13]. For e = 0, (16. 6ft gives an 
integral equation with a variational operator, which admits the same critical values and 
points as the differential one studied before. 

An efficient way to prove branching of solutions from p = n + 1, which applies for 
the lack of differentiability and regularity of nonlinearities involved, is using degree-index 
theory. This establishes branching from an isolated solution, say, the first one Fq for 
simplicity, from the branching point e = provided its index (rotation of the vector field 
/ - A'(F , 0)) satisfies [3H p. 353] 

(6.7) 7 = ind(F ,J-A'(F ,0))^0. 

Unfortunately, for the equation (16. 6p . this approach hardly applies, since the operator 
contains nonlinearities \F\~ a F, which are not differentiable at 0, so it is not completely 
clear how to treat the spectral properties of the self-adjoint operator G'(F ) for setting 
in the whole space E (not on a bounded interval corresponding to Dirichlet problems). 
Besides, the Frechet differentiability of A in F for e ^ is non-existent. Note that 
standard alternatives of bifurcation-branching theory without differentiability hypotheses 
assume sufficient regularity of the perturbations; cf. [TTj Thm. 28.1]. Justification of 
branching phenomena in the present problem needs further deeper analysis and more 
tricky "functional topology" involved. 

Continuing to describe spectral features of G'(F ), we note that 

A = (with the eigenfunction ip ~ Fq) 

is an eigenvalue that corresponds to the translational invariance of the original PDE (11.81) 
with the infinitesimal generator D y . Therefore, A = 1 is an eigenvalue of A'(F ,0) (this 
derivative exists), and this leads the critical case, where computing of the index is more 
difficult and is performed as in [311 § 24]. It is more important that 

A = 1 (with the eigenfunction ip\ ~ Fq) 

is also an eigenvalue of G'(F ) and this is associated with the generator D t of the group 
of translations in t. 
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Thus, the index condition (16. 7p needs special additional treatment, in particular, as- 
sociated with spectral properties of the linearized operator G'(F), and this is an open 
problem. In this connection, we conjecture that e = is a point of changing index for 
the basic family {Fi}, and that [31] Thm. 56.2] applies to generate a countable set of 
continuous e-curves from any basic pattern Fi(y) constructed in Section [531 

Continuing branching approach, for non-differentiable nonlinearities as in (16.6j) . as 
an alternative (and less effective relative to branches detected) approach, we have that 
Schauder's Theorem can be applied to get solutions of (16.51) for small e > and to trace 
out p-branches of the suitable profiles. Let us be more precise about this extension of 
p-branches from the variational critical points at p = n + 1. Thus, we consider (16.61) as an 
integral equation in R using the well-known spectral properties of the operator B* in Sec- 
tion [2] with compact resolvent; see [291 [22] f° r similar reductions. We then need to apply 
Schauder's fixed point Theorem [2, p. 90] in the framework of a weighted L^*-metric, in 
which the integral Hammerstein-type operators involved are naturally compact; see [311 
§ 17]. One can see that the right-hand side in (I6.6P is continuous in this topology at e — 0, 
so that this gives at least one solution which is close to Fq for e ~ 0. 

This somehow settles existence of at least one solution of (16 .4p for small |e| > in a 
convex neighbourhood of solutions Fi for p = n + 1 (i.e., e = 0). In this framework, 
uniqueness becomes a very difficult problem since the integral operators are not contrac- 
tive in this metric. But we then obtain at least a single continuous p;-branch emerging 
from p = n + 1. The global behaviour of p-branches is a hard problem to be tackled next. 

6.3. Bifurcation from constant equilibria. Here we study other bifurcation phenom- 
ena in this problem: the p-bifurcation from the constant equilibrium F(y) = for the 
ODE (16.41) (or (16.61) . where difficulties concerning non-differentiable nonlinearities do not 
occur. This leads to other non-basic profiles F. To this end, we use the known spectral 
properties of the adjoint operator (12.71) . Set 



where G(-) is the nonlinear part of the operator to be treated as a perturbation in the 



(6.8) F(y) = F* + Y(y), where F, = f? 




space L 2 p ,. We next set 
(6.10) 

Then (16.91) reduces to 



y 



az, where a A (3{l — a)\F* 



—at 



1 

r 



(6.11) 
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Y -a 4 G(z,Y,Y')- 



We finally apply classic bifurcation-branching theory [311 H5] for the equivalent integral 
operator with compact operators in L^»(R) constructed similar to (16.61) ; see also extra 
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details in [28l |9| [22]. It follows from the linearized operator in (16.111) that bifurcation in 
p can occur at the following points: 

(6-12) F(^I)ea(B*) = {-!} p = P« = n + l-f (1 > f). 

Recall that all the eigenvalues A; = — | of B* are simple, and hence, under typical assump- 
tions, correspond to bifurcation points; see e.g., [JX P- 381]. The rigorous justification 
in the framework of functional setting in weighted L 2 p * -spaces and, especially, asymptotic 
expansions as e = p — p\ — > become rather tricky and involved; we do not do this here 
and we refer to typical examples in [HI § 5]. This creates a countable set of p-bifurcation 
branches to be described next. The adjoint polynomials i[>*(y) do not change sign (see 
(12. 8p ). so that the patterns obtained for p rs p~t are expected to have a minimal number 
(sometimes, none) of intersections with the constant equilibrium F*. We do not know for 
sure how to classify these patterns and expect that as p — > n + 1 they are converted into 
the non-basic profiles F + 2k from Section [5791 with large k = k(l) > - (an open problem). 

Thus, in general, there exists a countable family of p-branches, and some of them, 
being extended to p = 1~, n = 0, are expected, after necessary scaling, to match with 
the countable set of linear profiles given in (|4.8p ; cf. [291 § 6.1]. A reliable rigorous 
identification of these p-branches for the integral equation of the type (I6.6P is a difficult 
analytical, as well as numerical problem. 

Since (I6.12p makes no sense for I = (and other small Z's), classic bifurcation theory 
[TT| p. 401] suggests that the basic p -branch of F (y) of the simplest shape (as well as 
other first ones with I < -) is assumed to exist for all p > 1. In other words, this branch 
cannot appear at bifurcation points such as those indicated in (I6.12p . In Figure [TBI we 
show the first p-branch of F (a) and the deformation of the profiles F (y), (b). We expect 
that this po branch is composed from stable solutions and hence represents the generic 
asymptotic blow-up behaviour for the parabolic PDE (jl.8p . 

A small part of the next px-branch of dipole-like profiles F\(y) for n = 1 andp e [2, 2.11] 
is shown in Figure [T71(a). where (b) demonstrates the corresponding deformation with p 
of Fi(y). Further extension of this branch beyond p = 2.17 leads to strong instabilities 
where the profiles suddenly jump to different shapes (which possibly belong to other p- 
branches nearby having the geometric structure as in Figure |9l we did not construct such 
neighbouring branches). 

The p-branches can connect various profiles, with rather obscure understanding of pos- 
sible geometry of such branches and their saddle-node bifurcations (turning) points. For 
p — n+ 1, the questions on connections with respect to regularization parameters are ad- 
dressed in [2"3"| § 7] posing problems of homotopy classification of patterns in variational 
problems and approximate "Sturm index" of solutions. 

For instance, Figure [TBI for n = 1 shows the connection of the profile F\(y) for p = 2.2 
(cf. the previous figure) with two "almost independent" mutually shifted profiles ±F for 
p = 2. Next Figure [191 in the enlarged form, explains formation of the zero set of profiles 
and shows in (b) that the eventual structure for p = 2 actually belongs for a member of the 
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family {— F , F }-gluing, i.e., it is the profile F_ 2 j j+2 (y), with exactly seven transversal 
zeros between ±F (y) structures. 

The principle fact that higher-order p-branches of the basic family {Fi} (see Section 
15.5)) can be originated from p = n + 1 is illustrated in Figure [2D] where we show the 
P2-branch of profiles F 2 (y) (a) (for p = n + 1, this profiles is given in Figure [7(c)) and the 
deformation of the profiles (b). 

Another principal issue is justified in Figure [21] that shows the p-branch of the profile 
^+2,2,4-2 (y); cf. Figure [8] for p = n + 1. 
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Deformation of F (y) into two almost independent F (y); n=1 




y 

FIGURE 18. Deformation of Fi(y) for p G [2, 2.2]; n= 1. 



Deformation of F (y) into two F (y): enlarged, n=1 _ 6 Deformation of F (y) into two F (y): enlarged, n=1 

D x 1 D 




FIGURE 19. Enlarged zero set for profiles in Figure [TBI (b) shows that for p = 2, 
it is F +2j 7-2(y). 



Concerning other, more complicated profiles for p — n + 1, such as F +i {y) and others 
containing such structures shown in Figures [S] and [1_21 numerical results suggest that these 
cannot be extended for p > n + 1. For instance, Figure 1221 demonstrates that the profile 
F +i {y) very quickly jumps to the type F +2 ,2,+2 (cf. Figure [21!) for the increment Ap = 
1CU 3 , and even for smaller Ap's. This illustrates the fact that both profiles are very close 
and are originated at the same branching point. We do not intend here to get numerically 
a correct bifurcation diagram for this delicate case. It seems that this demands more 

31 




p-branch of F for n=1 



Smooth deformation of F (y) for ps [2,2.05]; n=1 
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(b) F+2,2,+2 profiles 



FIGURE 21. The p-branch of solutions F+2,2,+2 of equation (|6.4p for n = 1 (a); 
corresponding deformation of -F+2.2,+2 (b). 

advanced parameter continuation techniques that are available in the standard MatLab 
environment. See a similar unstable case below for p < n + 1. 

Remark: on ^-bifurcations. The origin of various branches of solutions can be also 
seen via an additional parameterization; cf. [9j § 4.3]. Namely, we consider equation (16.9p . 
where we introduce the parameter by replacing 

(3 i-> fi > 0. 

Then, linearizing as above, we conclude that bifurcation points are (for even profiles) 
(6.13) 



1 

l ■ 



1 = 2,4, 
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Non-extendability of F +4 (y); n=1 
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FIGURE 22. F + 4 for p = n + 1 jumps to -F+2,2, +2 if is increased by 10~ 3 ; n = 1. 

Since /3 < 7 < |, we have that at least two bifurcation branches originated at /i 2 and 
/i4, we have a good chance to be extended to /i = (3 (this global continuation is an open 
mathematical problem, but can be checked numerically), and hence generate solutions of 
the original ODE fTO]) . 

7. Global blow-up similarity profiles for p g (l,n + 1) 

For p 6 (l,n + 1), we have /? < in (14. ip . so that the similarity solutions describe 
expanding as t — > T~ waves with global blow-up uniformly on any compact subset in x. 

7.1. Same oscillatory behaviour close to interfaces. The ODE (14.21) reads for / rs 

-(|/|»/) (4) -W + ... = 09 <0), 
so, replacing y — y f— > y, on integration for y w 0, we have 

(i/r/r=+/9w/+»-- 

This gives (13.31) . where A = /%o < is reduced to —1 by scaling. Thus, for p £ (1, n + 1), 
the similarity profiles are oscillatory near interfaces as for p — n + 1. By (13. 7p . the 2D 
asymptotic bundle is enough to match two symmetry boundary conditions (14. 4p . and the 
proof of matching remains open. 

7.2. Profiles and branches. It turns out that, for 1 < p < n + 1, the ODE (16.41) is more 
difficult to solve numerically than for p > n+ 1. In Figure [221 f° r n — 1, we present the 
Po-branch of the generic profiles Fq by continuation in p from p = 2.5 (single point blow- 
up) until p = 1.7<n + l = 2 (global blow-up). The actual continuous p-deformation 
of functions Fo(y) is seen from Figure [2H For smaller p > 1, the similarity profiles 
F (y) become larger (and according to (14.31) the corresponding constant equilibria diverge 
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p Q -branch for n=1 
5 1 1 1 1 1 ) 
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FIGURE 24. Deformation of F (y) from Figure [231 

exponentially as p — > l - ) and more oscillatory close to the interface. Note that the po- 
branch is expected to consist of asymptotically (structurally) stable blow-up profiles F (y), 
but we cannot prove this even in the linearized approximation (the linearized operator is 
a difficult non-self-adjoint operator with unknown spectrum). 

In order to avoid the exponential discrepancy (14.31) of branches as p — > 1~, we now 
perform in the ODE ( 16.41) the following additional scaling: 

(7.1) F^CF, ?/h ay, (7(i-°0(p-i) = a 4 = ^-(^^(p-i) (C = Fj^) 
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FIGURE 26. The p-branch of solutions Fz(y) of equation (|7,2p for n = 0.5 (a); 
corresponding deformation of -F+2.2,+2 (b). 



to get the equation with fixed equilibria F# = ±1 and 0: 

(7.2) -F® -(3(p-l)(l-a)\F\- a F'y-\F\- a F+\F\ p ( 1 - a ^ 1 F = («=^). 

In Figure l25| we show the po-branch for the ODE ( 17.21) with n = 0.5 and p G [1.1, 1.5]. 
In Figure |2"B1 we show the p 2 -branch of the profiles F 2 (y) = F„ 2 ,i, +2, 1,-2(2/) for the ODE 
(17. 2p with n = 0.5 (a small "discontinuity" is seen there). 

Concerning more complicated profiles not from the basic family {F/}, in Figures [271 and 
[28l we show the p-branches of the profiles F +4 (?/) and F +2i 2,+2(y) for the ODE (17. 2p with 
n = 0.5. We observe that both p-branches (and the corresponding deformations) look 

35 




very similar, and actually, they do coincide for 

p <p* « 1.487. 

We then claim that there exists a branching point p = p* G (l,n + 1), at which the p- 
branch splits into two, thus giving us two different profiles with similar geometric shapes; 
see [45] for classic branching theory. Then this is not a standard saddle-node bifurcation, 
and analytical (and numerical) justification of such branching for the equivalent integral 
equation with non-differentiable nonlinearities represents a difficult open problem. Note 
that, according to ( 16.121) . higher-order profiles F + 2k and geometrically neighbouring ones 
F+2,2,...,+2 can also appear at bifurcations from the constant equilibrium F*. This needs 
special research and remains an open problem. 

Finally, we must admit that our numerics did not detect the p-extensions to the global 
blow-up range p < n + 1 of the odd (anti-symmetric) basic profiles such as F\(y) ((b) 
in Figure [7j), or F 3 (y) ((d) in Figure [7j). We always observed a strong divergence of 
the method applied to the ODE (17. 2p without obvious reasons. Nevertheless, we expect 
that these p-branches of basic blow-up profiles do exist. This formal conclusion is associ- 
ated with the question of evolution completeness of blow-up solutions for the PDE (II. 8ft : 
indeed, then which patterns will describe a generic evolution as t — > T~ for classes of anti- 
symmetric initial data? Of course, these could be non-similarity solutions, but those hard 
questions are definitely out of the present rather preliminary self-similar consideration. 

Acknowledgement. The author would like to thank both the anonymous Reviewers 
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